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Abstract 

We look for differential equations satisfied by the generalized Jacobi polynomials 
|pQ,/3,A/,A'j-^^|^^^ which are orthogonal on the interval [—1, f] with respect to the weight 
function 

where a > -1, /? > -f , Af > and > 0. 

In the special case that f3 = a and N — M we find all differential equations of the 
form 



E 

i=0 



c,(x)y«(x) = 0, y(x)=P,r '''''' (^), 



where the coefficients {ci{x)}°^^ are independent of the degree n. 

We show that if M > only for nonnegative integer values of a there exists exactly 
one differential equation which is of finite order 2a + 4. 

By using quadratic transformations we also obtain differential equations for the poly- 
nomials \Pn'^^'^'^{x)\ for all a > —1 and > 0. 
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1 Introduction 



In the late thirties (see [|10|] and H.L. Krall classified all sets of orthogonal polynomials 
{P„(x)}^Q with degree[P„(x)] = n which satisfy a fourth order differential equation of the 
form 

Pi{x)y^'^\x) +P3{x)y^^\x) +P2{x)y"{x) +pi{x)y'{x) +po{x)y{x) = 

where {pi{x)}f^Q are polynomials with degree < i and {pi{x)}'^^^ are independent of 
the degree n. These sets of orthogonal polynomials include the classical Legendre, Laguerre, 
Hermite, Bessel and Jacobi polynomials. He also found three other sets of orthogonal poly- 
nomials satisfying a fourth order differential equation of this type. In Q A.M. Krall studied 
these new sets of orthogonal polynomials in more details and named them the Legendre type, 
Laguerre type and Jacobi type polynomials. These polynomials are generalizations of the clas- 
sical Legendre, Laguerre (with a = 0) and Jacobi polynomials (with /? = 0) in the sense that 
the weight function for these orthogonal polynomials consists of the classical weight function 
together with a Dirac delta function at the endpoint(s) of the interval of orthogonality. 

Later L.L. Littlejohn (see fl^) studied a generalization of the Legendre type polynomials 
and named them after H.L. Krall : the Krall polynomials. These Krall polynomials are 
orthogonal on the interval [—1, 1] with respect to the weight function 

^5{x + + ^6{x -1) + C, A> 0, B > and C > 0. 

In general {A B), these polynomials do not fit in the class of polynomials which satisfy a 
fourth order differential equation of the above type. The Krall polynomials satisfy a sixth 
order differential equation of a similar form. 

A.M. Krall and L.L. Littlejohn did some work on the classification of higher order differ- 
ential equations having orthogonal polynomial solutions. They tried to classify all differential 
equations of the form 

r 

5]k(x)2/«(x) = 0, re {2,3,4,...}, 

i=0 

where {Piix)}^^^ are polynomials with degree[pi{x)] < i and {Pi{x)}l^i are independent of n 
having orthogonal polynomial solutions {-Pn(a;)}^o with degree[P„(3;)] = n. See |14] and 

In T.H. Koornwinder found a general class of orthogonal polynomials which generalize 
the Legendre type, Jacobi type and Krall polynomials. These polynomials are orthogonal on 
the interval [—1, 1] with respect to the weight function 

r{a + P + 2) ^^^^^ ^ ^ ^^^^ ^ ^ ^^^^ _ 



2"+/3+ir(a + l)r(/3 + 1) 

where «>— 1,/3>— 1,M>0 and > 0. For these generalized Jacobi polynomials we will 
use Koornwinder 's notation : \ P"'^'*^'^(x) > 

L J n=0 ^ 

As a limit case he found the polynomials \ L^'^^ {x)> which are orthogonal on the 

I J n=0 

interval [0, oo) with respect to the weight function 

1 

x^e-"" + M6{x), a > -1 and M > 0. 



r(a + l) 

These polynomials generalize the classical Laguerre polynomials. 
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In P] J. Koekoek and R. Koekoek showed that the polynomials \ L'^'^ {x)\ satisfy a 

L J n=0 

unique differential equation of the form 

oo 

M^ai{x)y^'\x) + xy"{x) + (a + 1 - x)y'{x) + ny{x) = 0, 

i=0 

where {ai(x)}^Q are continuous functions on the real line and {ai{x)}°l^ are independent of 
n. It turns out that the coefficients {aj(x)}^Q are polynomials and the differential equation 
is of infinite order in general if M > 0. However, only for nonnegative integer values of a the 
order reduces to 2a + 4. 

We note that it is well-known that all sets of polynomials named before satisfy a second 
order differential equation with polynomial coefficients depending on n, but of bounded degree. 
See for instance 0], [13| and [|l5|] . 

In this paper we look for differential equations for the polynomials \ P"'l^'^'^ {x) \ with 
/? = a and N = M . Until now, only two special cases, due to H.L. Krall and L.L. Littlejohn, 
are known. In [|ic| ] H.L. Krall showed that the polynomials satisfy the 

following fourth order (if M > 0) differential equation : 

-^M(l - x^fy^^\x) + 4Mx(l - x'^)y^^\x) + 

+ (6M + 1)(1 - x^)y"{x) - 2xy'{x) + ^n(n + 1) [(n - l)(n + 2)M + 2] y{x) = 



and later L.L. Littlejohn found the following sixth order (if M > 0) differential equation for 
the polynomials |p^'^'^^'*^(x)| 



oo 
'n=0 



-^M(l - x2)3y(6)(^) _ Imx{1 - x^fy^^\x) - ^M(l - x^){l - 3x^)y^^\x) + 

+ ^Mx(l - x^)y^^\x) + (lOM + 1)(1 - x2)y"(2;) + 

- 4xy'(3;) + —n(n + 3) [(n - l)(n + l)(n + 2)(n + 4)M + 36] y(x) = 
36 

both in a different notation. The latter sixth order differential equation has not appeared in 
the literature yet. 

We will derive all differential equations for the polynomials \ P"'"'*^'*^(x) > for every 
a > — 1 and M > 0, which are of the form 

oo 

5:Q(x)y»(x)=0, 
1=0 

where the coefficients {cj(x)}^Q are continuous functions on the real line and {ci{x)}'^i are 
independent of the degree n. 

So, we consider the polynomials \ P"'"'^^'^^(x) > which can be defined by (see P], in a 

L J 71=0 

slightly different notation) 

pMM^^^ = CoP^'^)(x) - Cix^Pi°'-)(x), (1) 
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where 



Co = 1 + M- 



2n 



a + l) 



n 



+ 4M 



2/n + 2a + 1 



n 



(2) 



Ci 



2M /n + 2a^ 
{2a + 1) I n , 



+ 



2M2 



n + 2a \ I n + 2a + 1 



n 



1 



n 



(a + l) ^ 

As Koornwinder already remarked (see the case 2a + 1 = must be understood by 
continuity in a. 

Further we will show that for M > these differential equations are of infinite order in 
general and only for nonnegative integer values of a we find exactly one differential equation 
of finite order 2a + 4. This answers one of the questions raised in by W.N. Everitt and 
L.L. Littlejohn. 

Finally, we will also derive differential equations for the polynomials < '^^'^'^(a^) \ 
for all a > -1 and > 0. 



2 Some classical formulas 

In this section we give the definition and some properties of the classical ultraspherical poly- 
nomials |P7l°'"^(x)| ^. For details the reader is referred to Q and fl^ . We will only give 
those properties we need in this paper. Further we list some other classical formulas which 
we will use later on. ^ 

The polynomials \P^''^\x)\ can be defined by their representation as a hypergeo- 

L J n=0 

metric function as 

' 1 -x' 



n 



2 



, n = 0,l,2,.... (3) 



a + l 

A simple consequence of this definition is the differentiation formula 

L'*Pi"'")(x) 

/n + a\ / ly y^ (-n)fc(n + 2a + l)fc / l - x \^-' 
~ \ n ) \ 2) {k-i)\{a + l)u \ 2 ) 

- (";i(4y£^=%l^is^^(^)'-».^'^.---- 

Further we have the well-known symmetry relation 

P^--")(-a;) = (-l)"p(°.")(x), n = 0,1,2,.... (5) 

The ultraspherical polynomials satisfy a second order linear differential equation given by 

(1 - x^)y"{x) - 2{a + l)xy'(x) + n(n + 2a + l)y{x) = 0. 

By using induction it is easy to show that this differential equation implies that 

(1 - x2)D*+2p^°'")(x) - 2(a + i + l)2;D^+^P^"'")(a;) + 

+ (n-i)(n + 2a + i + l)D'P^"'")(x) = 0, i = 0, 1, 2, . . . . (6) 
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We also need the following formulas for the even and odd order ultraspherical polynomials : 



1\" /2n + Q^ 



n 



2F1 



-n, n + Q + 
1 
2 



x'' , n = 0,1,2, 



and 



respectively. These formulas can be found in |^ in a slightly different notation. 
We will often use the Vandermonde summation formula : 



2F1 



-n, b 
c 



1 



(C - b)n 



n = 0,1,2,... 



and the Saalschiitz summation formula : 



3-^2 



— n, o, b 
c, —n + a + 6 — c+1 



(c)n 

_ (c - a)n{c - b)n 
(c)n(c -a-b)n 



, n = 0,l,2,.... 



(7) 



, n = 0,l,2,..., (8) 



(9) 



(10) 



These summation formulas can be found in Q for instance. 

Finally, we remark that the Taylor series at the point zero of a hypergeometric function 
of the form p+iFp has a radius of convergence 1 unless it terminates. Moreover, such a series 
also converges absolutely at 1 if the sum of the numerator parameters is less than the sum of 
the denominator parameters. For details the reader is referred to ||l|. 



3 The differential equations for P° 



00 J 



In order to find all differential equations of the form 

00 

a,(x)y(*)(x) + (1 - x'^)y"{x) - 2{a + l)xy'(x) + n{n + 2q + l)y{x) = (11) 



i=0 



for the polynomials |P"'"'^^'^^(x)| ^, where the coefficients {aj(x)}^Q are continuous func- 
tions on the real line and where {ai(x)}^j^ are independent of n, we set y{x) = P"'"'^'^'^^(x) 
in ([Tl|) and use the definition @) to find 

00 00 
MCoJ2ai{x)D'P^'^^'^\x) - MC7i^ia^(x)L»^P^"'")(x) + 



i=0 



1=0 



MCixY,ai{x)D'+^Pj{''"\x) + 



1=0 



+ (l-x 

-2(q + 1)x 
+ n{n + 2a + 1 



+ 



+ 



C7oP^"'")(x) - Cix-^Pt'"\x) 



0. 
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Now we use (^) with i = and i = 1 to obtain 

OO r CO oo 

1=0 li=0 i=0 

We consider both sides as polynomials in M. Comparing the coefficients of equal powers of 
M on both sides leads to 



4 /n + 2a\^ 



M2 



(a + 1) 



n 



/ i=0 

n + 2a 



{2a + 1) \ n 



.1=0 



1=0 



M3 : 4 



4 /n + 2a\ /n + 2Q + lW2 
(oTT) \ n- 1 )\ n 

/ i=0 



n 
2 



n + 2a\ I n + 2a + 1 



(a + 1) I n - 1 



n 



.1=0 i=0 

It is clear that the third relation is a linear combination of the first two relations which can 
be simplified to 



5^ai(x)Z)*P^"'")(x) 



n 



1=0 



{2a + 1) I n I dx"^ 



and 



^ iai (x)Z^^Pi-'°) (x) + X a. (x)D'+ipi°'-) (x) = 4 



1=0 



1=0 



^n + 2a + l\ (f 
n — 1 / dx'^ 



(12) 



Pi"'")(x). (13) 



Since we demand that the coefficients {ai(x)}^^ are independent of n we introduce the 
following notation 

ao(x) := ao{n, a, x), n = 0, 1, 2, . . . 

(14) 

ai{x) := ai{a, x), i = 1, 2, 3, 

In order to find the general form of the coefficients {oj(x)}^Q we will prove the following 
theorems. 
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Theorem 1. The polynomials satisfy the infinite order differential 

L J n=0 

equation given by 



and 



^&,(x)y»(x) = 

1=0 

bo{x) := boin, a,x) = ^[l- (-1)"] , n = 0, 1, 2, 
2«— 1 

bi{x) := bi{a, x) = —-{-x)\ i = 1, 2, 3, 



(15) 



(16) 



Note that the coefficients {bi{x)}°^Q do not depend on a. The proof of theorem 1 can be 
found in the next section. 



Theore 



by 



m 2. The polynomials |p"'"'^'*^(x)| satisfy the differential equation given 

L J n=0 



Ci(x)y(*)(x) + (1 - x'^)y"{x) - 2(a + l)xy'{x) + n{n + 2a + l)y{x) = (17) 

i=0 



where the coefficients {cj(x)}^Q are defined by 

, , , , , , /n + 2a + 2\ 

co{x) := co[n,a,x) = 4{2a + 3)\ j , n = 0,1,2, .. . 



and 
where 



Ci{x) = {2a + 3)(1 - x2)c*(x), i = 1, 2,3, . . . , 
cl{x) := cl{a, x) = 



■■= c*{a,x) = ^J2 



2*^ / a + 1 \ A-2a-5\/l-x^'= 



A:=0 \ y 



A: y V 2 

i = 2,3,4, ... . 



(18) 
(19) 

(20) 



The proof of theorem 2 will be given in section 5. 

Now we will show that the general solution {ai{x)}°^Q of ( |T2|) and (13), where {aj(x)}^Q 
are continuous functions on the real line and where {ai{x)}'^i are independent of n, is 



ao(n, a, x) = ao{l,a, x)bo{n, a, x) + co(n, a,x), n = 0, 1, 2, 
ai{a, x) = ao(l, a, x)bi{a, x) + Ci{a, x), i = 1,2,3, ... , 



(21) 



where oo(l, a, x) is an arbitrary continuous function on the real line and where the coefficients 
{bi{x)}Zo and {ci{x)}Zo are given by (|T|), {^) and {^). 

The proof is based on the following lemma. 
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Lemma. Let {ai(x)}^Q be a set of continuous functions of the form which satisfies 
the homogeneous system 

( oo 

i=0 

(22) 

oo oo 
. i=0 i=0 

with ao(l, a, a;) = for a]7 reaJ x. Then we have 

ai{x) = 0, i = 0,1,2,... 

for a7] real x. 



In order to prove this lemma we substitute small values of n in the homogeneous system 
(|2^). Then we obtain for n = and n = 1 : 

ao(0, a, x)=0 and {a + l)[xao{l,a,x) + ai{a,x)] = 

for all real x. Since ao{l,a,x) = and a > —1, we conclude that ai{a,x) = for all real x. 
If we substitute n = 2 in (|2^) we obtain 



ao(2,a,a;)pf +a2(a,a;)^Pi°'"^(x) = 
2a2(a,x)^Pi"'")(x) +xao(2,a,x)^Pi"'"^(x) = 0. 



Hence 



ao(2, a, x) =0 



a2(a,x)^^P2"'"^(x) = -ao(2,a,x)P2"'"'*(x). 



Since 



1 



2Pi'^'")(x) - x^Pt'''\x) = -^{a + 2) ^ 



we conclude that ao(2,a,x) = for all real x and therefore that a2(a,x) = for all real x, 
since 



^i"'"^(^) = ^(« + 2)(2a + 3)^0. 



1 



dx^'' 2 
In the same way we obtain for n = 3 : 



3Pi"'")(x)-x^pi"'")(x) 



ao(3, a,x) =0 



a3(a,x)-^P3^"'"''(x) = -ao(3,a,x)P3^"'"''(x). 
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Now we find ao(3, a,x) =0 for all real x except for x = being the only zero of 
3Pt''"\x) - x-^Pt'"\x) = -i(a + 2)(a + 3)x. 

The continuity of oo(3,a,x) implies that ao{3,a,x) = for all real x. Then we also have 
03(0, x) = for all real x, since 

^Pt"'\^) = ^(« + 2)(a + 3)(2a + 5) / 0. 

If we proceed in this way we also find for each n > 4 that ao(n, a, x) = for all real x except 
for the possible zeros of nP^'°^^ (x) — x^P^"'"^ {x). The continuity of aQ{n, a, x) then implies 
that ao(^i a, x) =0 for all real x and finally we have a„(a, x) = for all real x, since 



dx" ^ ' \ n 

This completes the proof of the lemma. 

For the moment we introduce the term continuous sequence. A sequence {o-i{x)}^Q is 
called a continuous sequence if {aj(x)}^Q are continuous functions on the real line of the 
form (14). 

In section 4 (proof of theorem 1) we show that the continuous sequence {6i(x)}^Q defined 
by ([1^ ) is a solution of the homogeneous system (22). This implies that the general solution 
{aj(x)}^Q of (|2^), being a continuous sequence, is 

ai{x) = ao(l,a,x)6i(x), i = 0,1,2, 

This can be shown by using the lemma as follows. We define 

a*(x) := ai{x) — ao(l, a, x)6j(x), i = 0, 1, 2, . . . . 

Since ao(l, a, x) is continuous, {a*(2;)}^Q is a continuous sequence which is a solution of ( p2|) 
too where 

ao(l, a, x) = ao(li 2;) — ao(l, a, x)6o(l5 a, =0 

for all real x. Now the lemma gives us the desired result. 

In section 5 (proof of theorem 2) we show that the continuous sequence {ci(x)}^g given 
by dH), (ID and (H) is a solution of and (||). In order to prove that the general solution 
{ai(x)}^Q of ( [I2D and (p^), being a continuous sequence, is of the form (pi]), we simply note 
that {aj(x) — Cj(x)}^Q is a continuous sequence and a solution of the homogeneous system 
(ID. Hence 

ai{x) - Ci{x) = ao(l,a,x)6j(x), z = 0, 1,2, 

This shows that we have found all differential equations of the form (|^), where the 
coefficients {aj(x)}^Q are continuous functions on the real line and where {aj(x)}^^ are 
independent of n. 
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In M we conjectured that the polynomials |P"'"'*^'*^(x)| satisfy a differential equa- 
tion of the form (^) where the coefficients {cj(x)}^Q are given by (^) and (19) where 



4(-l)^/a + l\ / + + | 



x'\, i = 1,2,3, 



cl{a, x) = 



(23) 



-1)^+1 / a 



(2i + 1)! U - 1 



I (a + l)x 2F1 



-i + l,a + ^ — i 
3 ^ 
2 



i = l,2,3,, 

By using the formulas (0) and (^) we can prove that ( 22 ) is equivalent to 

c^(a, x) = 

c:(a,x) = ^p(_V+^'--+3)(x), z = 2,3,4, ... . 



(24) 



Here we remark that the parameter a — i + 3 in the Jacobi polynomial might be smaller than 
— 1. However, the ultraspherical polynomial Prl"'"^(x) is also a polynomial in a. Instead of 
the definition (^) we may define 



Pt'^Hx) = -L ^ ( (n + 2a + lUa + k + 



nl 



k=0 



X — 1 



n = 0,1,2, 



By using this definition of the ultraspherical polynomial we see that (24) is equivalent to (^ 
Since, by using (P), (|l|), ([l|) and (p|) we have 



(0) := a2i{a, 0) = (2a + 3) 



4(-l)^+^ (a + r 

m U-1 



y^O, i = 1,2,3,, 



if a is not a nonnegative integer we conclude that if M > the order of the differential 
equation given by ([Tl| ) can only be finite for nonnegative integer values of a. From the form 
( |23| ) we easily see that for nonnegative integer values of a we have 

degree[cj(2;)] = i, i = 2, 3, 4, . . . , 2a + 4 
Ci{x) =0, i > 2a + 4. 



Hence, for M > the order of the differential equation given by ( |llD can only be finite if we 
choose ao(l, a, x) =0 in view of (p]) and if a is a nonnegative integer. In that case the order 
of the differential equation given by (p^), (p^), ( p^ ) and (^) equals 2a + 4. 
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4 Proof of theorem 1. 

To prove theorem 1 we have to show that 

oo 



1=0 



i=0 1=0 

To do this we first note that we have in view of (^) 

oo -, oo 

Y^h{x)D^Pt''\x) = - [1 - (-l)-]P^")(x) +$:6.(x)D^p("'")(x) 

1=0 1=1 

and 

oo 1 I oo 

= - [1 - (-1)"] +^6,,(x)Z)^+ip^")(x). 



i=0 



i=l 



Now we use (I), (I) and dll) to obtain 



j=i 



1 /n + a\ ^ x^ (-n)fc(n + 2a + 1) ^ ^l-x^'^"^ 

1 /n + ^ (— n)fe(n + 2a + 1)^ ^ ^/c^ ^ ^ 1 — x 

2 i n 



E 

fc=i 
00 



E 

j=i 



1 / n + a \ ^ {-n)k{n + 2a + 1)^ 



2 \ n 



k=l 



{a + l)kk\ 



1 + X 



1 — X 



pf'K-x) - pf'Hx)] = "2 [1 - (-1)"] Pt'^'hx), 



Y^ih{x)D^P^^^^\x) 



i=0 



and 



I (n + a 
2\ n , 



E r-; _ m E 



{-n)k{n + 2a + l)fc / 1 - x 



(a + l)fc(A;-i)! 



k—i 



1 ( n + a\ ^ {—n)k{n + 2a + 1)^ ^ fk — l\ ^ fl — x 



E 



E 



H ^ y^tl (a + l),(fc-l)! f^,\i-lj 
l/n + a\ ^ (-n)fc(n + 2a + l)fc / 1 + x^^'^ 

(-l)"x-^Pi"'")(x), 
ax 



5:5,(x)Z)^+ip("'-)(x) 



i=l 



11 



1 / n + Q 



n 



1 / n + Q 



n 



1 / n + a 



n 



Sir i: 

j=l k=i+l 

(-re)fc(re + 2a + l)fc ^ 
^ (« + l),(A:-l)! ^ 

{-n)k{n + 2a + l)k 



-n)k[n + 2a + 1)^ / 1 — x 



(a + l)fc(/c-i-l) 



fc-i 



k=2 

oo 

E 

A;=2 



k-V 



X 



k-i-l 



1 — X 



k-i-l 



{a + l)k{k-l)\ 



l + x\^-^ fl-x 



This proves theorem 1. 



5 Proof of theorem 2. 

In this section we will give a proof of the main theorem 2. 
In view of (|l^) and (^) we have to show that 

and 



oo oo /_i_9_|_l\r/2 

j=0 j=0 — 1 



n 



Now we write by using (|lq), ([L9|) and (20) 



£ Q(x)Z)*Pi"'")(x) = 4(2a + 3) ( " + + ^ ) + £ Ci{x)D'P^^^^\x), 

CO oo 

Y^ic,{x)D^P^-'^\x) = Y,tc,{x)D^Pt''\x) 

1=0 i=2 

and 

oo / _1_ 9 _L 9\ ,^ °° 

y Ci(x)D^+ip("'")(x) = 4(2a + 3) r ^ -^P("'")(a;) + V Ci(x)D^+ip("'")( 

^ \^ n-2 ydx ^2 

Now we obtain by using (pC|) 

oo 

_ ^/g + A 2^+^ ^ (-2a + i-2)fc sfc p»+fc+2p(a.«). ^ 
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By using (^) we find 

(4) 



In + a 
I n , 



i+k+2 oo 



^ (-n)m(ra + 2a + l)m f l - x 



E 

m=k 



m - ky.{a + 1) 



m+j+2 



m—i—k—2 



m—k 



Now we use Vandermonde's summation formula (^) to see ttiat tliis leads to 



V (-2a + i - 2)k xfcp,i+fc+2p(a,Q)/ 



k=0 '"^ 

/n + a\ y^ y^ (-2a + i - 2)fc (-n)m+t+2(n + 2a + l)m+i+2 , . y+fc (1 
V " /toi; fcK^ + 3)fc (m-A:)!(a + lW,+2 2 



- x) 



2i+m+2 



' n + a 

\ ^ . 
' n + a 
I n 



' n + a 

' n + a 
I n 



(-2a + i - 2)fc (-n)m+i+2(ra + 2a + l)m+i+2 (-?7^)fc . . w.. „ 
fc!(i + 3)fc 2™+*+2(a + lW,+2 m! ^ 



m=0 fe=0 



m=0 



-n)m+i+2(?T- + 2a + l)m+^+2 (1 - X 

2i+™+2(ce + l)^+.+2 m! 



X 2i^l 



-m, —2a + i — 2 
i + 3 



.s^i >^ {-n)m+i+2{n + 2a + l)m+i+2 (1 - x)"^ (2a + 5)^ 
i^o 2™+*+2(a + lWi+2 m! (i + 3)^ 

^ j(i_+2)! ^ (-n)n^+i+2(ra + 2a + l)m+t+2 (2a + 5)^ 
" ) 2»+2 (m + i + 2)!(a + l)™-ui^9 m! 



'm+i+2 



1 — X 



m 



Hence 



^c*(x)Z)^p("'<^)(x) 



i=2 

oo 

E 

i=0 



' a 



+ 1^| - (-2« + ^-2). _ ..^m.+2p(„,„). 

(z + 2)!£:; fc!(. + 3). ^^"^ 



/n + a\ v^/_,y/ a + l\ (-n)^+■i+2(?^ + 2a + l)n^+^+2 (2a + 5)^ A - 
V /i^ V ("^ + * + 2)!(a + lW^+2 m! I 2 



' n + a 



Om+2("- + 2a + l)rra+2 (2a + 5)m ( \_ 

ml 



^ \-n)m+2[n + 2a + i)m+2 (2a + 5)™, / i - x 
("^ + 2)!(a + lW2 



X 



-n + m + 2, —a — l,n + 2a + m + 3 
m + 3, a + m + 3 



1 . 



Now we use the well-known fact that 



3-t'2 



a, 6, c 
d, e 



c — a 



-3^2 



a,b,c + 1 
d, e 



c — a 



-3F2 



a + l,b,c 
d, e 



z ] , c — a ^ 
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and the Saalschiitz summation formula dlG) to find for n > m + 2 



3F2 



-n + m + 2, —a — l,n + 2a + m + 3 
m + 3,a + m + 3 



1 



(n + 2a + m + 3) ^ f -n + m + 2, -a - l,n + 2a + m + 4 \ 
(2n + 2a + 1) ^ m + 3,a + m + 3 j"^ 

(n - m - 2) / -n + m + 3, -a - l,n + 2a + m + 3 
(2n + 2Q + 1)^ m + 3,a + m + 3 

(n + 2a + m + 3) (g + m + A)n-m-2{-n - 2a - l)n_m-2 _^ 
(2n + 2a + l) ("i + 3)„_m_.2(-"' - a)n-m-2 

^ {n-m-2) {a + m + 4:)n-m-3i-n - 2a)n-m-3 
{2n + 2a + 1) (m + 3)n-m-3{-n - a + l)n-m-3 
r(n + 2a + l) 1 (m + 2)! 



n! (a + m + 3) r(m + 2a + 4) 



[n(n + 2a + 1) + (a + l)(2a + m + 3)] . 



Note that the same result also holds for n = m + 2. For n < m + 2 we have (— n)m+2 = and 
the 3-^2(1) exists since 

(-n + m + 2) + (-a - 1) + (n + 2a + m + 3) < (m + 3) + (a + m + 3). 

This implies that 

00 

i=2 

_ fn + a] r(n + 2a + 1) -y^^ {-n)m+2in + 2a + 1)^+2 (2a + 5)^ 



n J nl m!(a + l)m+3 r(m + 2a + 4) 

X [n{n + 2a + 1) + (a + l)(2a + m + 3)] 



1 fn + a\ fn + 2a\ ^ {-n)m+2{n + 2a + l)m+2, , 

} ^ , (m + 2a + 4)x 

/ m=0 



(2a + l)4y n jy n J m!(a + 1)^+3 

X [n(n + 2a + l) + (a + l)(m + 2a + 3)] ( j . 

Since we have 

n(n + 2a + 1) + (a + l)(m + 2a + 3) = (a + m + 3)(m + 2a + 3) + (n - m - 2)(n + 2a + m + 3) 
and 

(m + 2a + 3)(m + 2a + 4) = m{m - 1) + 2(2a + 4)m + (2a + 3)(2a + 4) 

we obtain 

00 

1=2 

_ 1 fn + a\fn + 2a\ 

(2a + 1)4 V n )[ n ) ^ 



14 



E 



.m=0 



{-n)m+iin + 2a + 1)^+4 fl-x 
m\{a + l)m+4 



m+2 



+ 



+ 2(2a + 4)5] 



-"-)m+3("' + 2a + l)m+3 f - x\ 



m=0 



m!(a + l)m+3 



m+l 



+ 



+ (2a + 3)(2a + 4) > — 

^ ' ^„ m!(a + 1) 



)m+2 



+ 



m=0 



-n)m+4{n + 2a + l)m.+4 fl-x 
m!(a + l)m+4 



m+l 



+ 



-n)m+3in + 2a + 1)^+3 fl-x 



-(2a + 4)y 

m!(a + lW3 

n + 2a\ 



(2a + 1)4 \ n J 
■1-xV 



^^,p(-)(x)-4(2a + 4) ^ 



P^")(x) + 



+ (2a + 3) (2a + 4)^Pi"'")(x) + 
1 — 3; \ d' 



n + 2a ' 



(2a + 1)4 V n 



-(1 - x2)^Pi"'")(x) + 2(2a + 4)x^p("'-)(x) + 

^2 



+ (2a + 3)(2a + 4)— P("'")(x) 



Now we use (p) for i = 2 to find 



yc*(x)D^p(°'") 



i=2 



n + 2a 



(2a + 1)4 \ n y 

2(a + l)x^P^"'") (x) + [{n - 2)(n + 2a + 3) + (2a + 3)(2a + 4)] ^P^".")(x) 



Hence, by using ([1^) and (y) for i = 1 and i = we find 
J2ciix)D'Pt'''\x) 



i=2 



(a + l)(a + 2)(2a + l) 



(1 - X^) X 
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2{a + + [(n - 2)(n + 2a + 3) + (2a + 3) (2a + 4)] ^P^^^'^H 



(a + l)(a + 2)(2a + l) 



^n + 2a 
n 



4(a + l)(a + 2)x2-f-P^"'")(x) - (n - l)(n + 2a + 2)(1 - x^)— P^"'") (x) + 
- n(n - l)(n + 2a + l)(n + 2a + 2)P^"'")(x) + 



+ [(n - 2)(n + 2a + 3) + (2a + 3) (2a + 4)] (1 - x^)— (xl 



4 /n + 2a\ 
(2a + n Jdx^ 



.P^")(x) + 



n(n - l)(n + 2a + l)(n + 2a + 2) /n + 2a\ . ^ 

(a + l)(a + 2)(2a + l) i n I " ' 



(2a + 1) 
Finally, this implies that 



n + 2a + 2^ 
n - 2 , 



p(°'-)(x). 



oo 



^c,(x)Z)^p(°'-)(. 

1=0 



4 /n + 2a\ 



(2a + 1) I n 



which proves (25). 



The proof of (26) is much easier. We start from 

oo 

Y,tcl{x)D'Pt°'\x) 



1=2 

oo 

E 

1=0 



'a + l\ 2^+^ ^ (-2a + i-2) 



i,_n '*'-V' "T ^jfc 



We use Vandermonde's summation formula (^) again to obtain 

^ (-2a+J^-2)fc _ i+fc+2p(a,a)/ X 

_ /n + a\ (-2a + i-2)k {-n)m+i+2{n + 2a + l)n^+^+2 {-m)k , . 

[ n ) ^,to ^K^ + 2)fe 2-+*+2(a + lW,+2 m\ ^ ' 

i-n)m+i+2in + 2a + l)m+i+2 (1 - x)"^ 



'{l-x)"" 



n + a 

. n 



■m=0 



2'^+»+2(a + 1) 



'm+i+2 



ml 



X 2F1 



-m, —2a + i — 2 
i + 2 



' n + a 

' n + a 
I n 



>^ (-n)^+,+2(n + 2a + (1 - x)'" (2a + 4), 

- ) 2^ 2-+*+2(a + lWi+2 ml {i + 2)m 



m=0 



^w (i + l)! y> (-n)m+»+2(n + 2a + 1)^+^+2 (2a + 4)^ / I -x V 
" ^ 2^+2 t^ (m + i + l)!(a + l)„+,+2 m! I 2 J 
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This leads to 

oo 



i=2 

oo 



n + a\ a + l\ {-n)m+i+2{n + 2a + l)n^+^+2 (2a + 4 )^ fl-x 

+ a\ ^ (-n)m+2(n + 2a + 1)^+2 (2a + 4)^ fl-x 
\ ^ I ^0 {m + l)\{a + l)m+2 ml 



m=0 

X 3^2 



-n + m + 2, —a — l,n + 2a + m + 3 
m + 2,a + m + 3 



Now we use the Saalschiitz summation formula (|T0|) again to find for n > m + 2 

„ ( — n + m + 2, — a — 1, n + 2a + m + 3 
^ m + 2,a + m + ^ 

_ (g + m + 3)n-n^-2(-n - 2a - l)n-m-2 
(m + 2)n-m-2{-n - a)n-m-2 

r(n + a + 1) (m + 1)! r(n + 2a + 2) r(a + m + 3) _ (m + 1)! r(n + 2a + 2) 
~ r(a + m + 3) r(n) r(m + 2a + 4) r(n + a + 1) ~ r(n) r(m + 2a + 4) 

Since {—n)m+2 = and the 3i^2(l) exists for n < m + 2, this implies, by using (^) 

oo 
i=2 

fn + a\ ^ {-n)m+2(n + 2a + l)m+2 (2a + 4)„ (m + 1)! r(n + 2a + 2) fl-x 



^7,^0 ("^ + l)Ka + lW2 r(n) r(m + 2a + 4)V 2 



+ a\ r(n + 2a + 2) ^ {-n)m+2{n + 2a + l)m+2 fl-x 
^ n ) r(n)r(2a + 4) m!(a + 1)^+2 

4 /n + 2a + l\ d2 



(2a + 3) \^ n - 1 J dx"^ 
and therefore we have 

f:^c.{x)D^Pi-'-\x) = 4(" + + ^) (1 - x^)^p(-")(x). 

1=2 \ / 

Now we look at the second sum on the left-hand side of (p6|). Now we have 



Y^cl{x)D^^^Pt''\x) 

1=2 

^ /a + l\ 2^+' (-2a + i-2)fc _ ^ ,+^+3 



17 



As before we find by using tlie Vandermonde summation formula (| 



' n + a 

V ^ , 
'n + a 
i n 



^ ^ {-2a + i-2)k{-n)m+i+3{n + 2a + l)m+i+3{-m)k , .y+i/-, y 
to to + 2™W(« + lW,+3 m! ^ ^ ^ 

IC-lV+i V (-"-)rn+»+3(^ + 2a + l)m+i+3 (1 " a:)"" 



m=0 



m, —2a + i — 2 
i + 3 



/n + a\ ^ (-n)m+i+3(?^ + 2q + l)m+i+3 (1 - x)'^ (2a + 5)^ 

\ n P ' to 2-+i+3(a + l)„+,+3 ml (i + 3)^ 

(i + 2)! ~ ( 

1^ n ' 2^+3 ^0 (m + i + 2)!(a + l)™+i+3 ^^^! T 



■ X 



and therefore we obtain 

oo 

X:c."(i-)0'+'''^'"(a:) 



i=2 

1 fn + a 

2 i n 



X 



(-n)m+t+3(ra + 2a + l)m+»+3 (2a + 5)^ / I - x \"^ 
t {m + i + 2y.{a + l)rn+i+3 rnl \ o I 



1 /n + a\ ^ (-n)n^+3(?7- + 2a + 1)^+3 (2a + 5)^ / I - a^ V 
"2!^ n (m + 2)!(a + lW3 m! 12/ 



^ / — n + m + 3, — a — 1, n + 2a + m + 4 \ 

By using the Saalschiitz summation formula ( p!o| ) we have for n > m + 3 

p I — n + m + 3, — a — 1, n + 2a + m + 4 ^\ 
^ m + 3,a + m + 4: J 

_ (a + m + 4:)n-m-3{-n -2a- l)n-m-3 

(m + 3)„_m_3(-n - a)„_m-3 
r(n + a + 1) (m + 2)! T{n + 2a + 2) r(a + m + 4) _ (m + 2)! r(n + 2a + 2) 
~ r(a + m + 4) r(n) r(m + 2a + 5) r(n + a + 1) ~ r(n) r(m + 2a + 5) 

which gives us by using (^), since {—n)m+3 = and the 3^2(1) exists for n < m + 3 : 
00 



1=2 

_l/n + a\ {-n)m+3{n + 2a + l)nt+3 (2a + 5)r 
™ ) to {m + 2)\{a + l)m+3 ml 



X 
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(m + 2)! r(n + 2a + 2) fl-x 
r(n) r(m + 2a + 5) 



1 /n + a\ r(n + 2a + 2) (-n)m+3(n + 2a + l)m+3 / 1 - x 



2 1 n /r(n)r(2a + 5) 



E 

rra=0 



m!(a + 1) 



m+3 



(2a + 3)(2a + 4) 
This implies 



4 /n + 2a + l' 
(2a + 4) I n-1 ^ 



(2a + 4) ' ^ 
Now we have found 



n + 2a + 1' 
n — 1 



2(a + 2)x-^Pi"'")(x) - (n - l)(n + 2a + 2)^Pt°'\x) 



j=0 



j=0 



n + 2a + l' 
n-1 



(l-x2)^^'^"^(^) + 4(2a + 3) 



n-2 



dx 



+ 



4x 



(2a + 4) 



n + 2a + l 
n-1 



2(a + 2)x-^Pi"'")(x) - (n - l)(n + 2a + 2)-^Pi"'")(a 
ax^ ax 



n + 2a + 1 \ 



n — 1 



P(-'")(X) 



which proves (pq) and therefore theorem 2. 



6 Differential equations for P, 



a,±i 0,iV 



In this section we will derive a differential equation for the polynomials ^ Pn' ^'^'^ {x) 



and another one for the polynomials \ Pn '^ ''''^(x) 



n=0 



for all a > —1 and > 0. These 



n=0 



differential equations can be obtained from our results by applying the following quadratic 
transformations (see 0) : 



i-2r'^'^(^) 



Pn 



0,2M 



(2x2 - 1) 



and 



P 



a,a,M,M 
2n+l 



(x) xPn 



(2x2 _ 1) 



a,i 0,(4Q+6)Af 



pa,a,M,Mf']\ 



(27) 



(28) 
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Note that (|2|) and (|2|) reduce to (0) and (D if M = 0. 

If we set y{x) := f{2x'^ - 1) we can prove by induction that 

y'''\-) = E Tn^^^^.-''-''f'''(^-' - 1), . = 0, 1, 2, . . . 



{2j-2i)\{2i-j)V 



and 



2'+^ (2i + l)!23j-2*-i 

(2j-2i-l)!(2i-j + l)! 



or written in one formula 



i A\n3j-i 

y^'^(^)= E 7^7— TJ7— Tjx2^-7(^)(2x2-l), i = 0,1,2,.... (29) 



Note that (^) and (|29|) substituted in the differential equation given in theorem 1 leads 
to a triviality, since 

= y2^^-'x^^f^^\2x'-i)y- — ^-41 — - = o, 

in view of 

y LIT -i^y('^'V-iy-o ,-123 

^(2,-^)!(^-,)!- j! W ^ ^ -0,^-1,2,3,.... 
By using and ( p9[ ) we obtain from theorem 2 the following equation : 

+ (1 - [l6x2/"(2x2 - 1) + 4/'(2x2 - 1)] + 

- 8(a + l)xV(22;^ - 1) + 2n(2n + 2a + l)/(2x2 - 1) = 

with 

, , , , , , , /2n + 2a + 2\ 

co(x) := co(2n,a,x) = 4(2a + 3)l 2n - 2 I ' = 0' 1' • • • > 

satisfied by 

fix) := Pr'^''^''\x). 

Since 



20 



we obtain 



MY,dj{x)f^^\2x'^ - 1) + (1 -x^) \l6x^f"{2x^ - l)+4f'{2x^ - 1] 

j=0 

- 8(q + l)xV(2x2 - 1) + 2n(2n + 2a + l)fi2x'^ - 1) = 0, 



+ 



(30) 



where 



do{x) = 4(2a + 3) 



2n + 2a + 2^ 
2n - 2 , 



, n = 0,1,2,, 



2i 



{2j - - j)\ 



x^^-'ci{x), j = 1,2,3,.... 



By using ( p!9|) and (|2^) we easily see that dj{x) is an even polynomial with degree < 2j 
for each j = 1, 2, 3, . . .. Now we set 2x2 - 1 = t and = 2M in (|0|) to find a differential 
equation of the form 



NY,d*{t)y'^'\t) + {l-t')y"{t) + 

j=0 

- ^ [{2a + 1) + (2a + 3)t] y'(t) + ^n(2n + 2a + l)y(t) = 0, 

for the polynomials I Pn' ^'°'^(t) I , where 
I J n=0 



(iS(t) = TT(2a + 3) 



2n + 2a + 2' 
. 2n-2 , 



n = 0,1,2, 



and 



d*{t) = \dA J^^], j = 1,2,3, 



We remark that d*j{t) is a polynomial in t with degree d*{t) < j for every j = 1, 2, 3, . . .. By 
using (19) and ( p3D we see that 



d*i-l) 



]dm = ^^^^c,,{0) 



(2a + 3)(-l)^+^ 



'a + l\ 2^ 



7^0, i = 1,2,3, 



if a is not a nonnegative integer. This implies that the order of the differential equation is 
infinite in that case if > 0. For nonnegative integer values of a we have dj{x) = for 
j >2a + A and 

d2„+4(x) = 24«+8x2«+4c2„+4(x) / 

since degree[c2o+4(x)] = 2a + 4. This implies that the order of this differential equation equals 
2a + 4 if a is a nonnegative integer and > 0. 
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And if we set y{x) := xf{2x — 1) we find by using Leibniz' rule 



3=1 



(2j-2i)!(2i-j) 



2j-l 



(2i- l)!23i-2*+i 



(2j-2i + l)!(2i-l-j)! 



2i 

E 



(2i + l)!23i- 



■2j 



{2j-2i + l)\{2i-j)V 



and 



2j+l 

• E 



(2i + 1)!23J- 



2j-l 



(2j-2i-l)!(2i-j + l)! 



/^^^(2x^-l), f = 0,1,2,... 



+ (2i + 1) .„ :r^^-^7^^'^(2x^ - 1) 



(2j-2i)!(2i-j)!- 



(2i + 2)!23j 2» 1 (^.) 2 _ ^-012 

^(2j-2z)!(2z-, + l)!^ ^ 1), .-0,1,2,... 



or written in one formula 

(,) , ^ (i + l)!23^-- 



2/' 



= E 



{2j-i + l)\{i-j)V 



+'f(^\2x'-l),i = 0,1,2,.... 



(31) 



If we substitute (^ ) and (^) in the differential equation in theorem 1 we find a triviality. 



smce 



y(x) + E&,(x)y«(x) = xf{2x' - 1) + E E it ^ !^T~L ^''^'f'^'H'^^' - 1) 



i=l 



- (2j-i + l)!(i-j)! 

00 2j+] 

^23J-V^+VW (2x2-1) ^ 



(i + l)(-l)^ 



{2j -i + iy.{i-j)\ 



in view of 



2j+l 

E 



(i + i)(-i)^ 



i2j - i + iy.{i - jy. 



E 



;i + j + i)(-i)^+^' 



- (j-i + l)!i! 



-1)^ 



J+1 

E 

1=0 



3 

L 

i=0 



(-ir-E ' (-1) 



0, i = 1,2,3, 



By using (28) and ( |3lD we obtain from theorem 2 the following equation : 
Mf:c.{x) E .o/'T,Hwr .„ x^^-+V^^'^(2x2-l) + 



i=0 



(2j-i + l)!(i-j)!- 



+ (1 - x^) 16x^f"{2x^ - 1) + 12x/'(2x^ - 1) 



Vo™2 



+ 



2(a + l)x Ax^f{2x^ - 1) + 7(22;'' - 1) + 2(2n + l)(n + a + l)xf{2x^ - 1) = 
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with 

satisfied by 
Since 



cq{x) := co(2n + 1, a, x) = 4(2a + 3) 



2n + 2a + 3\ 
. 2n-l ) 



n = 0,1,2,... 



fix) := P:'^^'^^"'^'^'\x). 



i=0 



{2j-i + l)\{i-j)\ 



2j+l 



^23V(^) (2x2-1)^ 



(i + l)!2- 



j=0 

we obtain, after division by x 



{2j - i + - j)] 



x^^-'+^Ci{x), 



M^ej(x)/(^)(2x2 - 1) + (1 - [l6xV"(2x2 - 1) + 12/'(2x2 - 1) 

3=0 



+ 



- 2{a + 1) 4xV(2x2 - 1) + /(2x2 - 1) 
+ 2(2n + l)(n + a + l)f{2x'^ - 1) = 0, 



+ 



(32) 



where 



2j+l 



+ 1)!23J- 



; (2j-z + l)!(^-i)! 



x^^-'aix), j = 0,1,2,.... 



By using ( [T^ and ( p3[ ) we easily see that ej(x) is an even polynomial with degree[ej(x)] < 2j 
for each j = 0, 1, 2, . . . . Now we set 2x2 - 1 = t and iV = (4a + 6)M in (|3|) to find a differential 
equation of the form 

oo 

NYe*it)y^'Ht) + {l-t^)y"it) + 

j=0 

- i [(2a - 1) + (2a + 5)t] y\t) + ^n(2n + 2a + 3)y(t) = 0, 
for the polynomials < Pn ' ^ ' ' (i) f , where 



em 



n=0 

1 /2n + 2a + 3\ 

2 \^ 2n - 1 J 



, n = 0,l,2,... 



and 



em 



1 / ll + t 



J = 1,2,3, 



8(2a + 3)''^ 

Note that e*(t) is a polynomial in t with degree ej(i) < j for every j = 1, 2, 3, . . .. By using 
( p^ ) and(p3D we see that 

^a + 1^ 



e-i-l) 



1 . , 2^~^ 

e,(0) = (-1)^+^' 



8(2a + 3) 



' a 



(a + 1) + {2j + 1) I 



J-l, 



(2a + 5)(-iy+i ^ ^ / 0' J = 1, 2, 3, . . . 
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if a is not a nonnegative integer. This implies that the order of the differential equation is 
infinite in that case ii N > 0. For nonnegative integer values of a we have ej{x) = for 
j > 2a + 4 and 

e2a+4{x) = 24-+«x2"+^C2«+4(:r) / 0. 

This implies that the order of this differential equation equals 2a + 4 for nonnegative integer 
values of a and N > 0. 
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